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Domain walls formed by one dimensional array of vortex lines have been recently predicted to 
exist in disordered helical magnets and multiferroics. Using a long wavelength approach supported 
by numerical optimization we lay out detailed theory for dynamics and structure of such topological 
fluctuations at zero temperature in presence of weak disorder. We show the interaction between 
vortex lines is weak. Also we show the internal degree of freedom of this array leads to the enhance- 
ment of its mobility. We present estimates for the interaction and mobility enhancements using the 
microscopic parameters of the system. 

PACS numbers: 75.10.-b, 75.60, 75.70, 75.85 



I. INTRODUCTION 
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Domain walls and domain patterns in magnetic ma- 
terials to a large extent determine the controllability of 
the magnetic behavior of devices, a crucial factor in tech- 
nology of today's spintronics and magnetic storage sys- 
tems. Amongst various types of magnetic structures, he- 
lical magnets are in current interest demanding detailed 
studies of their properties under various circumstances. 
With the discovery of multiferroics 1 , materials such as 
RMn0 3 in which R G {Y^D^Ho} 2 .^ with coupled 
magnetic and electric properties and their intrinsic spiral 
magnetic order the necessity of such studies is under- 
standable more than ever. A new class of topological do- 
main walls has been predicted to exist in Helical magnets 
in recent studies^. These domain walls separate two anti- 
chiral domains and posses vorticity. These domains are 
indeed observed in circularly polarized X-ray scattering 
experiments^. Static domain patterns in these obser- form- 
vations suggest strong pinning of the walls. In another 
recent studyii the pinning mechanism of domain walls 
in these systems (Hubert walls) has been presented. Fol- 
lowing this study, with an emphasis on pure vortex walls 
we present a detailed theory for structure and dynam- 
ics of these walls. We show that inside the topologically 
protected vortex walls^ the interaction between vortices 
is weak. We explain the effect of disorder on their struc- 
ture and dynamics. We also explain how the internal 
dynamics of such walls can affect their mobility in pres- 
ence of weak disorder. This is indeed a new realization of 
previously considered systems^ where internal degree of 
freedom enhances the mobility of the magnetic structure 
and therefore reduces the threshold force density. 



Wc consider a major class of helical magnets, the cen- 
trosymmetric (CS) systems in which the magnetic mo- 
ment distribution is invariant under space inversion. The 
helical axis in CS ground state is fixed. The Ginzburg- 
Landau hamiltonian density for helical CS magnets fre- 
quently used in literature can be written in the following 




FIG. 1: Schematic illustration of the cross section of a vortex 
line (oval region) formed between two anti-chiral domains. 
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h{m] = - (<9 2 m + q 2 m) + (V±m) 



(1) 



in which the helical axis is dictated by the RKKY ex- 
change anisotropy to be the x axis. 



II. VORTICES IN HELICAL MAGNETS 

The existence of vortex domain walls in Helical mag- 
nets and multiferroics indicates the importance of a de- 
tailed analysis of vortex fluctuations in such systems. 
These systems are highly anisotropic as the the system 
develops a ferromagnetic order in each easy plane and 
a helical order in the perpendicular direction. It is still 
possible to have vortex line excitations with the vorticity 
inside the easy planes. These vortex lines would consist 
of conventional vortices in each plane and their energy 
per unit length would be e± — J/a\ogL/a in which L 
is the size of the ferromagnetic planes. On the other 
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hand there are vortex line excitations in which the di- 
rection of the vortex line is parallel to the easy planes 
and the structure of these vortices is not as simple as 
the former ones. In this case one needs to solve the non- 
linear saddle point differential equation associated with 
the hamiltonian (pj along with the boundary condition 
V x Vip = 27r<5 2 (r)z. This has been performed before 7 
using an ansatz and separation of scales. In the next 
section we show that their result is also confirmed with 
numerical minimization. The variational approximation 
is used for different scales using a function of the form 
cp = tan (Ay/x) which minimizes the energy by varying 
A however we do this only for a certain distance r defined 
as: r 2 = x 2 + X(r) 2 y 2 . The result is 7 -: 
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e*(r) = 4 log 1 / 2 (r/a) 
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9 2 log(r/a) 



1/2 



(2) 



and: 



A 2 (r) 



5/64 
log(r/a) 



(3) 



For large enough r ^> ae ' 08 ^ 2 (the constant is not 
important since it comes from a variational calculation) 
the logarithmic term in the bracket dominates and the 
energy behaves logarithmic similar to conventional vor- 
tices in two dimensional XY models. For scales much less 
than this, close to the core of the vortex the value of the 
energy behaves as square root of logarithm. As the result 
this type of vortices in helical systems posses a layered 
structure. 



III. VORTEX WALL 



Another topological excitation possible in helical sys- 
tems is vortex domain wall which is a one dimensional 
lattice of vortex lines parallel to the easy plane axesi. 
These walls separate antichiral domains from each other 
(see figure Q] for better imagination) . Using the same 
ansatz before one can introduce the following ansatz for 
the vortex wall: 
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n=-N v 



tan 



Ay 



(4) 



in which x n = nir/q and 2N V + 1 3> I is the number of 
vortex lines. This solution indeed satisfies the boundary 
conditions tp(x, y) = ±qx and d y ip = as y — > ±oo. This 



can be seen by noting that for N v 3> 1: 
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(x - x n ) 2 + X 2 y 2 



N v 
n=-N v 

iq , , . 
= — (cot z — cot z) 

_ X{x — x n ) 

- {x _ Xn) 2 + A 2 y 2 

(cot z + cot z) 



n=-N v 

qX 
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(5) 



(6) 
(7) 



in which z = q(x + iXy). Clearly the above two expres- 
sions satisfy the boundary conditions. 

From the above calculation it is interesting to observe 
that the field distribution close to the vortex wall is iden- 
tical to the field of a single vortex: 



d x ip(z -> 0) 

dy<f(z -> 0) 
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x 2 + X 2 y 2 

Xx 
x 2 + X 2 y 2 



(8) 



indicating that the value of the variational parameter in 
smaller scales will be identical to the value obtained for 
a single vortex but this would not be the case for larger 
length scales (see below). Moreover this is indicating the 
fact that the individual vortices in a vortex wall almost do 
not disturb each other in this variational solution, hence 
their interaction must be weak. This will be confirmed 
also numerically in the next section. It is very instructive 
to derive an analytic expression for the deformation en- 
ergy of such vortex lattice as we proceed. The GL Hamil- 
tonian for the low energy fluctuations of the vortex wall 
will be explained by m(r) = ei cos ip(r) + e 2 sin <p(r) and 
their energy density, using ([1]) is: 
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<d 2 V ) 2 \ (9) 



In presence of the vortex wall at larger length scales we 
can write the effect of the wall deformations in the field 



as: 



d x ip = A x (r) +tp 



(10) 



in which A x (r) = qsgn q [y — Uy(x.\\)] is a smooth sign 
function which varies on the scales of q~ x and it y (x||) is 
the position of the wall at each point of its plane and 
ip is the remaining of the effects which is supposedly 
small. As the result up to second order in if) we can 
write [(d x ip) 2 — q 2 ] ~ Aq 2 ip 2 . However we re-insert the 
full ip field into the Hamiltonian in order to preserve the 
non-perturbative effect of the vortices in the field: 



{(d ± ipf+6 2 [d x <p-A x (r)] 2 y (11) 
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Here we have ignored the highest derivative of tp assuming 
d x ijj <C qijj. Introducing the rescaled coordinates f = 
{x/9,y, z} and fields A = 9{A X ,0, 0} we can write the 
above energy in the form of the vortex line lattice energy 
in a superconductor: 



J 



Vip - A(r) 



The saddle point equation from the above energy: 



V<p-A 



= 



must be solved in presence of the vortices: 
V x W(p = J(r) 

in which: 

J(r) = 2tt J2 / dr «' 53 ( r - r n) 

n=-N v •* 



(12) 



(13) 



(14) 



(15) 



The solution to the above system can be obtained by 
introducing a vector potential: 



Vif = A + V x a 



(16) 



with a gauge freedom which we use by choosing V • a = 0. 
Inserting this solution into (|14p we obtain the equation: 



V 2 a = V x A- J 
which has the solution: 



(17) 



a(r) =J G(r -r')[V' x A- J(f')] (18) 
Inserting this solution into the energy we obtain: 

' / [V x A - J(f )] • [V' x A - J(f ')] (19) 
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Let's calculate the vector A = V x A — J more explicitly. 
In order to do that we use the large scale approximation 
in which the density of the vortex lines is uniform : 

26» 3 

A(f) = -^-{d z u y y + i)[5 q {y - u y ) - 5{y - u y )] (20) 

this function is obviously non-zero only at a narrow re- 
gion with the thickness q^ 1 of the order of the core of the 
vortices where our continuum approximation fails. The 
equation (TT9|) clearly shows that the helical background 
represented by the vector potential screens the vortic- 
ity fields in large scales consistent with our choice of the 
ansatz ©. 



FIG. 2: Left: Vector field illustration of numerically opti- 
mized interface between two ani-chiral phases. Here 8 = 0.32. 
Right: Same as left but tilted. The plus/minus sign indicates 
the change in chirality hence existence of Hubert wall. 
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FIG. 3: Numerical comparison of the energy of the vor- 
tex domain wall per vortex(circles) with the results of the 
ansatz(solid line). The numerical result is obtained by opti- 
mization of the discrete model. 



IV. NUMERICAL INVESTIGATION 

For a more realistic investigation of the properties of 
vortex domain walls in helical magnets it is possible to 
use optimization techniques to find the lowest energy pos- 
sible configurations of the discrete model at zero temper- 
ature. It is also possible to investigate the interaction 
between vortices in such configurations. In this section 
we show that the ansatz solution for the vortices obtained 
previously are in good agreement with optimized results 
in this section. Also we show that the interaction be- 
tween vortices in a pure vortex wall are comparatively 
weak in agreement with our previous analytical specula- 
tions. We also justify the theory^ for the tilt of vortex 
domain walls via generation of staircases of Hubert walls 
and vortex lines. 

We use a two dimensional system with the assumption 
that the system is symmetric in the third direction i.e. 
the direction of vortex lines. We start with a configura- 
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FIG. 4: Vortex interaction energy as a function of distance 
from the wall. 



tion consisting of two helical domains with opposite chi- 
rality with a sharp interface between them and find the 
lowest energy possible using local iterative optimization 
methods. We use periodic boundary conditions in the 
direction of helical axis thus the number of the vortices 
in the wall is always even. For a wall with n vortices(see 
FigEJ we find the total energy per unit of the length 
of the wall, l n = nir/9 (in units of lattice constant and 
J = 1) obeys the following form: 



a + — A/log 4, 



(21) 



fitting gives j3 « 3.54. This gives the energy of one vortex 
to be proportional to (logn/O) 1 / 2 apart from an unim- 
portant constant in agreement with the ansatz (Figj3]). 

In order to make an estimate of the order of the in- 
teraction between these vortices in the wall we take the 
numerically obtained solution and displace one vortex 
out of the wall and obtain the minimum solution again 
by the same optimization technique(Fig|4])- It is inter- 
esting to see that the energy difference between the two 
configurations is almost vanishing. This is a strong evi- 
dence that the vortices in the vortex wall almost do not 
interact and without pinning centers to hold the walls 
the domains may not exist. 

The most elementary excitation of such wall based on 
the above results then are the formation of staircases 1 -. 
These staircases are consisting of vortices with a Hubert 
wall in between. It is easy to see this structure formation 
by considering a tilted array (see figure [5] right). 



TILTED DOMAIN WALLS: 
FORMATION 



STAIRCASE 



We showed in sections IIVI and IIIII that the vortex lines 
in a vortex wall interact very weakly. A tilt of the vor- 
tex wall must then generate a staircase of vortex lines 



of width ir/q. Recently a theory for elasticity of vor- 
tex walls has been developed^ based on vicinal surface 
approach. In this approach the elasticity of the vortex 
wall is found out to be local. Here for the sake of com- 
pleteness we explain in more detail the derivation: Let's 
assume the tilt is around an axis parallel to ferromagnetic 
planes. This means the generated vortex lines are par- 
allel to each other. The density of vortex lines for such 
a tilt with an angle a is n(a) = (g/7r)tana. It is also 
important to remember that the energy of vortex lines in 
the system is distributed highly anisotropic (see seclHl 
Eq. [2]) so the energy of each vortex line would depend on 
a. As the result by adding the energies for the Hubert 
segments and the vortex lines (neglecting the vortex in- 
teractions) the energy density per unit of the area of the 
wall will be: 



a(a) — 



e v (a)n(a) + (T H 



cosa| ( — e„(a)tana + an 



(22) 



(23) 



In which cry is the elastic constant of the Hubert wall 11 . 
Note that the energy density of a vortex line e v depends 
on a because of its anisotropic nature. A small tilt e away 
from this configuration will cost an extra energy density 
to create vortex lines: 



Aa(a, e) 



a(a + e)L(a + e) — a{a)L{a) 
a'(a)\e\ + ^[cr(a) + a"{a)}e 2 



(24) 



where L(a) is the linear size of the wall. The elastic 
energy of this wall then would be: 



n e (a) 



J d 2 r||Aa(a,e(r||)) 



(25) 



Note that the linear term is zero at finite temperatures 
because the roughening transition temperature is almost 
vanishing 11 . Above this temperature the vortices slide 
against each other almost freely so the average e will be 
zero. With this fact we conclude that the total energy of 
the staircase is: 



H e (a) 



\c\{a)(d x uf + c 2 (a)(d z ii) 



(26) 



which is local unlike Hubert walls 11 . Note that the term 
(d z u) 2 exist because of the area change of the wall upon 
tilting around x-axis. The elastic constants depend on 
the tilt of the domain wall and can be numerically calcu- 
lated. The results are shown in figure [5j 



VI. DISORDER 

We consider a major type of disorder in this paper 
namely defects in crystals where magnetic atoms are 



5 



cifa) 




jt/2 ji 



FIG. 5: Elastic constant of the vortex wall as the result of 
vicinal fluctuations as a function of orientation for 9 = 0.5. 

missing on a particular site or a non-magnetic impurity is 
replaced by the magnetic atom. These defects are called 
vacancies. Usually domain walls gain exchange energy 
by adopting to these sites. In this case we attempt to 
derive a Hamiltonian for these types of disorder in con- 
tinuum limit. Domain walls experience friction force due 
to these types of disorder. External driving forces can 
mobilize these walls provided they are stronger than a 
threshold value called depinning threshold. Periodic na- 
ture of the vortex walls however can have interesting be- 
havior in presence of disorder and driving force. The 
internal lattice movement of the wall helps avoiding the 
pinning centers hence increasing the mobility of the walls 
above threshold value. Threshold values for both vortex 
wall and Hubert wall have been presented elsewhere^. 
Here we explain the general theory in detail and show a 
new effect, the enhancement of the mobility of the wall 
solely due to its internal degree of freedom. 

A. Dynamics of discrete domain walls 

Here we discuss the dynamical aspects of such domain 
walls far from equilibrium in presence of external forces. 
We also show how it is possible to extract certain char- 
acteristics of domain walls such as Larkin length and de- 
pinning threshold. We show how a non-uniform density 
of a domain wall will result in the increase of effective 
mobility compared to generic values. Let's start with 
the overdamped equation of motion for each vortex line 
in the array: 

M- 1 -^ = -|^ + f(x,u„)+h (27) 
at ou„ 

in which u„ = {u n (z), v n (z)} is the displacement vector 
of the vortex line number n with y and x components 
and: 

f(x,u„) = -^ (28) 

is the friction force density as the result of disorder. Here 
x = {a;, z} representing a point on the wall. Also in the 
above Tl is the elastic Hamiltonian of the wall determined 



by equation (|26l) however throughout the rest of this cal- 
culation we assume the most general elastic Hamiltonian: 

n = Jj k \u k \ 2 (29) 

in which * indicates the Fourier transformation. Finally 
h is the deriving force density presumably applied uni- 
formly across the wall. To determine Hi, the disorder 
energy we use the fact that the wall consists of a one 
dimensional periodic array of vortex lines: 

Hi = Jk{t)£{t) (30) 

in which: 

£ (r) = ^2 ev[x-nn/q-v n (z),y-u n (z)] (31) 

n=-N v 

and e v (x, y) is the energy density of a vortex line at point 
r = {x,y, z}. The function JC(r) determines the effect of 
disorder: 

K(r) = 7 ^[/ ; (r-r,)- n/ ] (32) 

i=l 

with //(r) being the vacancy /impurity form factor 
J r fi( r ) — 1 an d nj is their density. Throughout this 
article we assume the form factor vanishes in a range 
of the order of the crystal lattice constant. The sum in 
the above is over the position of impurities/vacancies 
which is random, 7 is the effective volume occupied by 
impurity/ vacancy , Aj is the number of impurity sites 
and n./ = Ni/£l the impurity concentration. The above 
energy can be written as: 

Hi = fC(r)^2e v [x-rm/q-v n (z),y-u n (z)] = 

Jr n 

= / K,(r)e v (x - x',y - y')p(r') 

J r.r' 

= Jv(v) Pv (v). (33) 

where we have defined: 

V(v) = J^lC(r')e v (x-x',y-y')6(z-z') (34) 
p v (r) = S[y - u n (z)]8[x - rnr/q- v n (z)] (35) 

n 

as impurity potential and vortex line density. Here we 
approximate the vortex line density as if the displacement 
of vortex lines vary smoothly over the domain wall as the 
result we can treat the displacement vector as a smooth 
function of x: u n (z) ~ u(x, z). The impurity energy will 
be then of the following form: 

Hi = d 2 x V(x,u v )p(x — u x ) (36) 
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where we have defined p(x) = ^2 n S(x — rm/q). The 
force density vector has components perpendicular and 
parallel to the domain wall plane f = {/||, fx}' 

dV 

/j_(x,u) = -—p(x-u x ) (37) 

/||(x,u) = V(x,u v )p'(x-u x ) (38) 

The goal is to find the force density h c at which the wall 
starts to move. However application of the most conve- 
nient tool namely perturbation in powers of disorder will 
only give us the force-velocity relation of the system at 
high velocities. While this approach will not ultimately 
lead to a relation in the full range of external force mag- 
nitudes it provides us with an estimate for the effective 
mobility of the system and also it will ultimately help 
us to estimate the threshold force density in large length 
scales in a renormalization group point of view-i2. Let's 
now use the perturbation theory order by order to clarify 
the above argument: First we introduce the average drift 
velocity of the wall v which in our case has only a non- 
zero component perpendicular to the domain wall plane 
v = vy. The fluctuations out of the flat plane then can 
be separated: 



u(x,i) = vi + £(x,i) 



(39) 



On the other hand the disordered medium creates static 
corrugations with period <7 _1 in the wall: 



(€(x,t))=g(x) 



(40) 



Here: 



(Mr)) = / n 



tA'j 



i-O-^p.r.) 



means averaging over disorder. This effect does not exist 
in more conventional Neel or Bloch type domain walls but 
here the periodic nature of the wall dictates such changes. 
This effect however is conventional in studies of pinned 
lattices^. We will present the resulting corrugation later 
in this section. 

Now we expand the out of plane displacement in orders 
of disorder strength: 



(0) 



;(2) 



(42) 



after averaging for the zeroth and first order we will have: 



A 4 1 ' v 



and: 



|«(k, W ) = 
4 1} (k, W ) = 



£(0) = 



^V(k x + Q,k v ±)p{-Q) 
<> ipj 1 ^ + Sj^k) 
^V{kx + Q,ky,%)K-Q) 



(43) 
(44) 



(45) 



(46) 



ip^ 1 uj + d| 1 (k) 
in which we have assumed the elastic energy for fluc- 
tuations parallel and transverse to the domain wall are 
governed by two distinct functions G\\ and Q± in the form 
of equation |2"91 For second order on the other hand: 



f(x,u) =f(x,vf)+£- V u f(x,vt) 



(47) 



(41) and after averaging: 



p^v — h 



- J GZ\* - xOfff (x') + (fM ■ V u / X (x, vt) 



(48) 



the second term on the right hand side depends only on integrating over one period using the fact that g][' (0) 
x and is periodic consequently depends only on x 0: 
and is periodic as well. Taking Fourier transform and 



-(p/v - h) 



iQyQx 



<?x X ( k ) ~ ip^QyV Qn x (k) - ipu X Q V 



in the above we have used the correlator function: 

A(r) = (V(v)V(O)) - 7 2 n / /,(r) 



|p(Q x )| 2 A(k,Qj,) 



(49) 



(50) 
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where: 

/„(r) =5(z) J dx'dy'e v (x-x',y-y')e v (x',y') 
The depinning threshold naturally follows from setting v — > 0: 



k,Q 



iQ 3 y g±(k) + iQyQlGn (k) \p(Q x )\ 2 A(k x - Q x , k z , Q y ) 



g± (x) A< 3 ) (x, 0)C(x) + G» (x)A« (x, 0)D(x) 



(51) 

(52) 
(53) 



in which: 



(54) 
(55) 



and the derivatives of the correlator are introduced as 
follows: 



C{x) = / p(x')p(x - x') 

J x' 

D(x) = / p'(x')p'[x' -x) 



AW(i 
A^(r 



(d y V(r)V(0)) 
(d*V{r)d y V(0)) 



(56) 
(57) 



Here we try to explain in details an approximate method 
used before^ to estimate the threshold force density. The 
first assumption is that the correlator of disorder poten- 
tial A(r) and its spatial derivatives are of short range 
with an average value of Av' (for nth derivative) at 
a length scale L. The length scale dependence comes 
from a renormalization point of view where one system- 
atically sums up the contributions from smaller scales 
into the correlator using methods such as perturbation 
etc.— With this assumption and taking into account the 
fact that D(x), C(x) > we arrive at the estimate for 
the threshold of force density at the length scale L: 



G±(L)A® +G\\(L)A 



(58) 



On the other hand for an isotropically distributed set of 
pinning centers, A^ is generically zero for odd n re- 
sulting in the zero threshold force density however it is 
well knowr>i2 that in a renormalization group flow A^ 
develops a pole at scales larger than Larkin length scale 
C where the energy gain due to collective weak disorder 
fluctuations can overcome the elastic deformation energy 

( 2^ 

cost. This means A^ ' (x, it) at L w C has a cusp near 

{x, it} = resulting in A^ being nonzero. The second 
term from the parallel deformations however remains zero 
at all length scales. Using this relation then one can esti- 
mate the threshold force density without having to worry 
about the physical value of the mobility. This will give 
us the maximum value possible because for scales higher 
than the Larkin length the disorder energy gain tends to 
decrease although it stays higher than deformation en- 
ergy cost. 



The above argument also leads us to a way to deter- 
mine the Larkin length scale of the problem. Here we 
ignore the in-plane fluctuations of the wall and treat the 
vortex wall as a wall with uniform density. This will still 
give us the upper bound for the threshold force density 
because the internal fluctuations of the wall tends to re- 
duce the magnitude of the threshold force as we will see 
at the end of this section. We start with the one-loop 
renormalization group flow equation for the force corre- 
lator R(x,u) = (/x(x,u)/ x (O,0)) = -A( 2 )(x,w). Pro- 
vided we are close to the upper critical dimension we can 
use the renormalization group equation for Rl—' 



dR L (u) 
dL 

A(L)5L 



A(L)—[R L (u)(2R L (0) 

^ 2 (k) 



Rl(u))} 



(59) 



taking two times derivative of the first equation in (|59")l 

(3) 

and setting u = and assuming R L (0) =0 as an ana- 
lytic even function, we find the following solution up to 
the zeroth order in e: 



[#i 2) (0)] 



[4 2) (o)] 
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X ' L d?k 



l/a 



87T 



2 <? 2 (k) (60) 



in which R a (u) is at L = a the bare correlation function. 
As was mentioned earlier this solution is divergent at the 
Larkin length given by the following equation at zeroth 
order of e: 



i?i 2 )(0) = 



^S 2 (k) 
l/o 8?r 



-1 -1 



(61) 



In order to find the renormalized non-zero value of 
R'(0) = A( 3 )(0) at Larkin length scale we need to use 
the functional renormalization group equation again. As- 
suming A{L) = AqL 1 ^ 1 we use the following ansatz: 



R L (u) = aL^R+ifiuL-^) 



(62) 



The parameter a is fixed by the condition that the fixed 
point function R*(Q) = 1. Inserting this into the RG 
equation results in a = — 2£ + 7 and: 

1/2 

R' L (0) = L-^[R L (0)}^ ( ^P-) ' (63) 



2An 



8 



On the other hand using the fact that J_°° Ri,{x)dx is 
invariant under the flow we obtain the exponent £ with 
the physical assumption that R*(x) — >• as x — >• ±00 to 
beC = 7/3. 



B. New effects in mobility correction 

For the mobility itself on the other hand an interesting 
effect reveals when one calculates the effect of fluctua- 
tions of the wall internal structure. Using the real part 
of the equation l|49p one can find a correction to mobility 
at non-zero velocity: 



v = jl±h 



1 



(64) 
(65) 



in which: 



ID 



t± = 



Q v \p(Q x )\ 2 A(k x + Q x ,k z ,Q y ) 

(66) 



£± f Q 2 y Ql\p(Q x )\ 2 A(k x + Q x ,k z ,Q y ) 

Mil A,Q 



^(k) +{^ x vQ v f 



(67) 



Using the fact that A(k, Q) — r y 2 ni\i v (kx,Q)\ 2 (e v does 
not depend on k z ) is always positive one can see that 
£± < corresponding to the out of plane fluctuations of 
the wall trying to reduce the effective mobility of the wall 
while £ii > tries to increase it, an effect which depends 
on the gradient of the density of the wall (note the term 
\Qxf>(Qx)\ 2 in the second equation). Also £\\ = for 
uniform density. This is physically plausible since the 
lattice structure of the wall while driven, helps the wall 
to avoid resisting centers during its motion. Note that 
when fin — >• we have £n — > 0. 

We can estimate the value of the above corrections by 
using proper cutoffs for large distance approximations 
and at high enough velocities. Those considerations leads 
us to the following approximation in units of lattice con- 
stant (a = I): 



C. Corrugation of The Vortex Wall 

In this section we briefly analyze the structural prop- 
erties of the domain wall in presence of domain wall. 
Our starting point is the equation (|48l) from which we 
can deduce the average displacement of the vortex wall 
at higher velocities. Since the second term on the right 

(2} 

hand side of this equation is periodic the function g]_ (x) 
will have to be periodic as well. For simplicity we assume 
Qx = a±k 2 and Q\\ — cr\\k 2 . After this consideration and 

some algebra one obtains for g±(x) = ^2 m gme lmqx the 
following: 



^(1 



Jm.O) 



njtO 



III- 



q — T^d n 



(70) 



in which: 



-ik y A(k) 



& j_7t 2 7 k er±M± [k 2 + {mq - k x ) 2 ] - > ny 



ik„v 



dm — 



p\\q 



-ik y A(k) 



tt 2 ct±_ 7 k <J\\p\\ [k 2 + (mq - k x ) 2 } - ik y v 



(71) 



Interestingly we see that the above coefficients are all 
zero for v — which means the corrugations appear or 
enhance when the wall starts to move. It is also possible 
to estimate the value of the above integrals at the moving 
phase of the wall. Using e„(k) = —Jq/ak 2 we will have: 



c±i « d±i 



71/7 



1/2 



(72) 



where we have the same simplifications as previous sec- 
tion for mobility and elastic constants. Also in the above 
we have used a cutoff for the integral over k x to include 
only wavelength of the order of q^ 1 . Finally smaller 
wavelengths (\m\ > I) in the corrugation have been ne- 
glected for large scale approximations. 



VII. DISCUSSION AND CONCLUSIONS 



■II 



a n [ — ] q 

-nn [—) q 



(68) 
(69) 



In the above we have used Q x = mode of density to 
estimate the correction to mobility coming from trans- 
verse fluctuations and Q x = q to obtain the effect of the 
internal degree of freedom. Also we have made simplifi- 
cation by assuming p = fj,± « /in and used the previous 
approximations 7 a± ss mi « Jq/a in a quadratic elastic- 
ity form. 



In this paper we have presented an analytic model 
as well as numerical verification of the predicted vortex 
walls in both clean and disordered helical magnets. From 
an analytic point of view we showed how the helical back- 
ground screens the vorticity effect in long distances re- 
sulting in the weak interaction between vortex lines. We 
have also presented a detailed discussion of the dynamics 
of such wall in the presence of disorder. Using established 
rcnormalization group ideas we have explained in detail 
how it is possible to theoretically predict the threshold 
force density for domain walls in general. In this cal- 
culation we have also shown that a moving vortex wall 
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array will display on average spatial corrugations. In this 
discussion we have also revealed how the periodic na- 
ture of the wall will result in correction of the mobility 
toward higher values compared with other conventional 
magnetic domain boundaries. Here we can further dis- 
cuss such effect by time scale considerations. The time 
scale during which a domain wall with a thickness of 
~ q^ 1 passes a pinning center during its motion is of 
the order of t± ~ q~ 1 /v. On the other hand the effect 
of this impact throughout the wall propagates during a 
time scale T|| ~ \ 2 /ayi^ in which A is the length of prop- 
agation. The ratio of the two time scales then will be: 
T ||/ T -L ~ #(V£||) 2 m which £|| — {aw^a/v) 1 / 2 is the typ- 
ical correlation length between fluctuations throughout 
the wall. In order for the periodic nature of the wall 
have any effect in the dynamics of the wall the propa- 
gation length must be at least of the order of the pe- 
riod of the vortex line array A ~ q^ 1 . These consid- 
erations result in a range of velocities v < v p ~ crnfinq 



at which the periodic effect show up. This will then 
directly determine the range of external force densities 
h < hp = v p /fj.± ~ t7||<7/i||//i_i_ at which the periodic 
structure affects the motion of the wall. This value of h p 
can be used in experiments to study the nature of the 
boundaries between anti-chiral domains. 
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